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Abstract. In this paper we study the homogenization of unsteady 
Stokes type equations in the periodic setting. The usual Laplace op- 
erator involved in the classical Stokes equations is here replaced by a 
linear elliptic differential operator of divergence form with periodically 
oscillating coefficients. Our mean tool is the well known two-scale con- 
vergence method. 



1. Introduction 

Let Q be a smooth bounded open set in (the iV-numerical space M. N 
of variables x = (xi, ...,xjy)), where N is a given positive integer, and let T 
and e be real numbers with T > and < e < 1. We consider the partial 
differential operator 

^ dxi \ l °d Xj/ 

1,3=1 J 

in Ox]0,T[, where a\- (a?) = ay (f) (x £ fl), a tj 6 L°° (R^;R) (1 < i,j < 
N) with 

(1.1) dij = dji, 

and the assumption that there is a constant a > such that 

N 

(1.2) a v M OCi > « ICI 2 for all C = (C 3 ) € R N and 

for almost all y € R^, where R^ is the iV-numerical space R^ of variables 
y = (yi, ...,7/at), and where |-| denotes the Euclidean norm in R^. The opera- 
tor P £ acts on scalar functions, say ip £ L 2 (0,T; i? 1 (O)). However, we may 
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2 LAZARUS SIGNING 

as well view P e as acting on vector functions u = G L 2 (o, T; F 1 (fi) w 
in a diagonal way, i.e., 

(F E u) 8 = FV (i = l,..,iV). 

For any Roman character such as i, j (with 1 < i,j < N), u l (resp. 
u 3 ) denotes the i-th (resp. J-th) component of a vector function u in 

Lj" oc (flx]0,T[) JV or in L\ oc (R^ X R r ) where R T is the numerical space 
R of variables r. Further, for any real < e < 1, we define u £ as 

u £ (x,t) = u(j,^\ ((ar,t) G Ox]0,T[) 

for u G ^ oc ( M ^ x as is customary in homogenization theory. More 

generally, for u G L/ oc (Q x Ry x R T ) with Q = fix]0,T[, it is customary 
to put 

u e (x, t)=u fx, t, ~, ~j ((s, t) G Ox]0, T[) 

whenever the right-hand side makes sense (see, e.g., [SJ 9]). 

After these preliminaries, let f = (/*) G L? ($,T;H- 1 (fi;R)^. For any 
fixed < e < 1, we consider the initial boundary value problem 

(1.3) + P £ u £ + gradp E = f infix]0,T[, 

(1.4) drou e = in Ox]0,T[, 

(1.5) u e = on dQx]0,T[, 

(1.6) u £ (0) = in n. 

We will later see that as in [18J, (jl.3|) - (|1.6p uniquely define (u £ ,p £ ) with 
u e G W (0, T) and p e G L 2 (0, T; L 2 (ft; R) /R) , where 

W (0, T) = {u G L 2 (0, T; V) : u' G L 2 (0, T; V') } 

V being the space of functions u in i^o 1 (OjR)^ with alivu = (V is the 
topological dual of V) and where 

L 2 (O; R) /R=|i)g L 2 (O; R) : J 

Let us recall that W (0, T) is provided with the norm 

i 

ll U llw(0,T) = (|I U IIl2(0,T;U) + ll u/ !lL2(0,T;y')) " ( u ^ W (0, T)) , 

which makes it a Hilbert space with the following properties (see [18]): 
W (0, T) is continuously embedded in C ([0, T] ; L 2 (0) N ) and is compactly 

embedded in L 2 (o, T; L 2 (ft) 7 



valx = 

n 



UNSTEADY STOKES TYPE EQUATIONS 



3 



Our aim here is to investigate the asymptotic behavior, as e — > 0, of 
(u £ ,p e ) under the assumption that the functions a%j (1 < i, j < N) are peri- 
odic in the space variable y. The steady state version of this problem (i.e., 
the homogenization of stationary Stokes type equations) was first investi- 
gated by Bensoussan, Lions and Papanicalaou [2]. These authors use the 
well-known approach of asymptotic expansions combined with Tartar's en- 
ergy method. We also mention the paper of Nguetseng and the author [16] . 
on the sigma-convergence of stationary Navier-Stokes type equations. 

The present study deals with the periodic homogenization of an evolution 
problem for Stokes type equations. 

This study is motivated by the fact that the homogenization of (|1.3p - (jl.6p 
is connected with the modelling of heterogeneous fluid flows (see, e.g., |19| 
for more details about such models). 

Our approach is the well-known two-scale convergence method. 

Unless otherwise specified, vector spaces throughout are considered over 
the complex field, C, and scalar functions are assumed to take complex 
values. Let us recall some basic notation. If X and F denote a locally 
compact space and a Banach space, respectively, then we write C (X; F) 
for continuous mappings of X into F, and B (X; F) for those mappings 
in C (X; F) that are bounded. We shall assume B (X; F) to be equipped 
with the supremum norm \\u\\ = sup^g^ ||w(x)|| (||-|| denotes the norm 
in F). For shortness we will write C (X) = C (X; C) and B{X) = B (X; C). 
Likewise in the case when F = C, the usual spaces LP (X; F) and L p loc (X; F) 
(X provided with a positive Radon measure) will be denoted by LP (X) and 
L p loc (X), respectively. Finally, the numerical space M. N and its open sets are 
each provided with Lebesgue measure denoted by dx = dx\...dxN- 

The rest of the paper is organized as follows. Section 2 is devoted to the 
preliminary results on existence and uniqueness, and some estimates on the 
velocity u e , the pressure^ and the accelaration of the fluid, whereas in 
Section 3 one convergence theorem is established. 



Let f2 be a smooth bounded open set in R , let T > be a real number 



and let f = (f j ) € L 2 Jo.T;^ 1 (Sl) N ). For < e < 1, it is not apparent 



that the initial boundary value problem (jl.3p - (jl.6p has a solution (u e ,p £ ), 
and that the latter is unique. With a view to elucidating this, we introduce, 
for fixed < e < 1 the bilinear form a 6 on H% (Q; R) N x Hi (fi; R) N defined 
by 
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for u = (u k ) and v = (v k ) in (fi;R) . By virtue of (jl.ip . the form a E is 
symmetric. Further, in view of (jl.2|) . 

(2.1) a e (v,v) >a||v||^ 1(n)JV 

for every v = G i^Q 1 (ftjR)^ and < e < 1, where 




dx 



with Vi> = yjj^, On the other hand, it is clear that a constant 

Co > exists such that 

(2.2) \a e (u, v)| < c ||u|| H i (n) iv 

for every u, v G Hq (0; R)^ and < e < 1. 

We are now in a position to verify the following result. 



Proposition 2.1. Suppose f iies in L 2 ^0, T; L 2 (fi; R) j. Under the hy- 
potheses il.l\) - [TJ^) , the initial boundary value problem tl.3\ )- [lJ)\) deter- 
mines a unique pair (u e ,p e ) with 

u £ G L 2 (o, T; H% (O; R) n ) n C f[0, T] ; L 2 (O; R)^) and 



p e G L 2 (0,T;L 2 
Proof. For fixed < e < 1, we consider the Cauchy problem 

f u' £ (t) + A £ u £ (t)=£(t) in]0,T[ 
{2 - 6) U(0) = 0, 

where A £ is the linear operator of V into V defined by 
(A £ u, v) = a 6 (u, v) for all u, v G V 
and t is the function in I? (0, T; V') defined by 

(£ (t) , v) = (f (t) , v) for all v G V 

and for almost all t g]0, T[, and where (,) denotes the duality pairing be- 
tween V' and V as well as between H~ l (f2;R) and i?Q (fijlR)^. Thanks 
to (12.1j) - (12.2p the Cauchy problem (12.3P admits a unique solution u £ in 
>V(0,T), as is easily seen by following [H Chap. 3, Theoreme 1.2, p. 116], 
see also [18} pp. 254-260]. Now, let us check that the abstract parabolic 
problem fj2.3[) is equivalent to (|1.3p - (|1.6p . Let U e (t) = J Q P e u £ (s)ds and 
F (t) = Jq f (s) ds for < t < T, where u e satisfies (|2.3p . It is evident that 
U e and FeC ( [0, T] ; H' 1 (O; M)^) . By the first equality of ([13]) we have 



(2.4) j t (u £ (t) ,¥») = (£ (t) - A e u e (t) , ¥») for all ^ G V, 
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where 



V = |v? G T> ($7; R) JV : dwy> = j . 

Integrating (|2.4p . we have 

(u e (t) + U £ (t) - F (t) , if) = 0, < t < T, if G V. 

Thus, using a classical argument (see, e.g., [18j p. 14]), we get a function 
P £ E C ([0, T] ; L 2 (0; R) /M) such that 



Hence p e = G X>' (Q) and the pair (u £ ,p £ ) verifies fjl .3|) (in the distribu- 
tion sense on Q), with in addition (|1.4p - (jl.6p . of course. Futhermore, by us- 
ing the fact that f G L 2 (o, T; L 2 (O; R) N ^ we have u' £ G L 2 (o, T; L 2 (O; R)*) , 

as is easily seen by following [HI p. 268] . Therefore p £ lies in L 2 (0, T; 1? (Q; R) /R) 
and is unique. Conversely, it is an easy exercise to verify that if (u £ ,p £ ) is a 
solution of ([IU-([L6]) with (0, T) and p e G L 2 (0, T; L 2 (0; R)) , then 

u e satisfies ()2.3p . The proof is complete. □ 

The following regularity result is fundamental for the estimates of the 
solution (u £ ,p £ ) of (fL3]) - (fT6|) . 

Lemma 2.1. Suppose f, f" G L 2 (0, T; V) and f (0) G L 2 (Q; R)^. T/ien f/ie 
solution u e o/ | QL 3|) verifies: 



where H is the closure ofV in 1? (il;R) . 

The proof of the above lemma follows by the same line of argument as in 
the proof of |18|. p. 299, Theorem 3.5]. So we omit it. We are now able to 
prove the result on the estimates. 

Proposition 2.2. Under the hypotheses of Lemma 1 2. 11 there exists a con- 
stant c > (independent of e) such that the pair (u £ ,p £ ) solution of \1.3\j - 
fOI) in W (0, T) x L 2 (0, T; L 2 (O; R) /R) satisfies: 

(2.5) ll u ellw(0,T) — c 



u £ + U £ + gradP £ = F. 



G L 2 (0, T; V) n L°° (0, T; H) 



(2.6) 



< c 



dt 



L 2 (0,T;H-^{Q.) N ) 



and 
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for almost all t G [0, T]. By taking in particular v = u e (t) in (12. 8p . we have 
for almost all t G [0, T] 

j f |u e (i)| 2 + 2a K (i)f < ^ ||f + a K (i)|| 2 

where |-| and ||-|| are respectively the norms in I? (Q) N and H% (Q) N . Hence, 
for every s G [0, T] 

\u E (s)\ 2 + a [ S \\u £ (t)\\ 2 dt< - F \\f(t)\\ 2 v ,dt 
Jo a Jo 

since u e (0) = 0. By the preceding inequality, we see that 

(2.9) a [ T \\u £ (t)f dt < - [ T \\f (t)\\ 2 v , dt. 

Jo « 7o 

On the other hand, the abstract parabolic problem (]2.3p gives 

u ; = f-,4 £ u e . 

Hence, in view of (12. 2p 

(2-10) ll U £lL 2 (0,T;V") - H f Hl 2 (0,T;V') + C \\ u s\\ L 2 (0,T;V) • 

Thus, by (|2.9|) and (|2.10p one quickly arrives at ()2.5|) . Let us show (|2.6|) . 
We are allowed to differentiate (j2.8j) in distribution sense on ]0, T[, and by 
virtue of the hypotheses of Lemma |2. H we get u" G I? (0, T; 1/') and 

(2.11) «,v)+a £ (i4,v) = (f',v) (v G V) . 

In view of Lemma 12. 11 we take in particular v = (i) in (|2.8 j) . This yields 

|<(t)| 2 + a £ (u £ (t) X (*)) = (f (*)) (* 6 [0,21) • 

Further, since G L 2 (0, T; V) and u£ G L 2 (0,T;y'), w e have G 
C ([0, T] ; H). Hence, by taking in particular t = in the preceding equality 
and using (jl.6p one quickly arrives at 

X (0)| 2 <|f(0)|X (0)|, 

i.e., 

(2.12) X(0)|<|f(0)|. 

The inequality (|2.12p shows that u^. (0) lies in a bounded subset of H. On 
the other hand, by taking in particular v = u^. (t) in (|2.1ip . we get 

| X (t)\ 2 + 2a |X (t)\\ 2 < \ V (t)|£, + a |X (t)|| 2 

for almost all t G [0, T]. Integrating the preceding inequality on [0, t] (t G 
[0,T]) leads to 

ft 



X W| 2 + ot I |X (*)|| 2 ^ < - ll f 1lL 2( oT-yo + K (°) 
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It follows from (|2,12p and the preceding inequality that u' £ belongs to a 
bounded subset of L 2 (o, T; L 2 (O; M) N j . Hence f|2.6j) is immediate. Let us 

prove U22J. For almost all t G [0,T], p £ (t) G L 2 (fi;R) /M. Thus, by US p. 
30] there exists v e (i) G (fi; M) N such that 

(2.13) divv £ (t)=p £ (t) 
(2-14) ||v e (i)||< Cl |ife(*)l 

where the constant c\ depends solely on Q. Multiplying (II. 3p by v £ (t) yields 

« (i) , v E (i)) + a £ (u e (t) , v E (i)) - f p £ (t) divv e (t) dx = (f (t) , v £ (t)) 

is 

for almost all i G [0, T]. Integrating the preceding equality on [0, T] and 
using (|2TT3>(j2T4) lead to 

II^IIl2(Q) ^ ^HI^L^T;!/^^) H^Hl^Q) + C l ll f llL 2 (0,T;H-i(n)) llPe IL 2 (Q) 
+CiC ||u e || L 2( 0)T .y) ||P £ |Il 2 (Q) , 

where c is the constant in the Poincare inequality, cq and c\ are the constants 

in (|2~2|) and ([234]) respectively. Thus, 

(2-15) 

lbe||l,2(Q) < c l c || u e|L2( 0iT;i 2( C )iv^+Ci ||f|lL2(0,T;H- 1 (n))+ c l c II u £|Il 2 (0,T;V) • 

Combining (f2TT5l) . (|23|) and (gj)]) leads to ([277]) . □ 

3. A CONVERGENCE RESULT FOR (}1.3|) - (|1.6|) 

We set y = (— 5, f ) , y considered as a subset of (the space M. N 
of variables y = (y\, j/at))- We set also Z = (— 3,5), Z considered as 
a subset of M T (the space R of variables r). Our purpose is to study the 
homogenization of (jl.3jl - (jl.6p under the periodicity hypothesis on etjj. 

3.1. Preliminaries. Let us first recall that a function it G Lj oc (M.y x R T ) 
is said to be Y x Z-periodic if for each (k, I) G Z N x Z (Z denotes the 
integers), we have u(y + k, r + l) = u(y,r) almost everywhere (a.e.) in 
(y, r) G 1^ x R. If in addition u is continuous, then the preceding equality 
holds for every (y, r) G R N x R, of course. The space of all V x Z-periodic 
continuous complex functions on x M T is denoted by C per (Y X Z); that 
of all Y x Z-periodic functions in Lf oc (R^ x R T ) (1 < p < oo) is denoted by 
Lper (Y x Z). C per (X x Z) is a Banach space under the supremum norm on 
R N x R, whereas Lp er (y x Z) is a Banach space under the norm 

u(y,T)\ p dydr) P (u G L p per (Y X Z)) . 



Lp(VxZ) 
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We will need the space H\ (Y) of Y-periodic functions u G LL\ oc {^y) 
^loc fay) such that f y u (y) dy = 0. Provided with the gradient norm, 



\u\\ H i #{Y) = ^\V y u\ 2 dy^j (ueH^(Y)), 



where V y u = (jjyj, g|^J> H# (Y) is a Hilbert space. We will also need 
the space L^ er ( Z; (Y) J with the norm 



L$ er (Z;H^(Y)) = [ j v l Y \VyU{y,T)\ 2 dydT^ (u G L 2 per (Z '; H\ (Y))) 



which is a Hilbert space. 

Before we can recall the concept of two-scale convergence, let us introduce 
one further notation. The letter E throughout will denote a family of real 
numbers < e < 1 admitting as an accumulation point. For example, 
E may be the whole interval (0,1); E may also be an ordinary sequence 
(£n) ne N w ith < e n < 1 and e n — > as n — > oo. In the latter case E will be 
referred to fundamental sequence. 

Let be a bounded open set in R% and Q = Qx]0,T[ with T G W + , and 
let 1 < p < oo. 

Definition 3.1. A sequence (u £ ) eeE C L p {Q) is said to: 

(i) weakly two-scale converge in L p (Q) to some uq G L p (Q; L v ver (Y x Z)) 
if as 

E3e^0, 
(3.1) 



u e (x, t) ip e (x, t) dxdt uq(x, t, y, r) ij) (x, t, y, r) dxdtdydr 

J J JQxYxZ 

for all V G L p ' (Q; C per (Y X Z)) (jr = 1 - ±) , where ^ (x, t) = 

^(M,f,£) ((x,t)eQ); 

(ii) strongly two-scale converge in L p (Q) to some uq G LP (Q; Lp er (Y x Z)) 
if the following property is verified: 

Given n > and v G L p (Q; C per (Y x Z)) with 
\ U Q ~ v \\lp(QxYxZ) — 2' ^ iere * s some a > such 
that \\u £ — v £ \\ LP ^Q-j < 7] provided E B e < a. 

We will briefly express weak and strong two-scale convergence by writing 
u £ — > uo in L p (Q)-weak 2-s and u £ — > uq in LP (Q)-strong 2-s, respectively. 

Remark 3.1. Lt is of interest to know that if u £ —¥ uq in L p (Q)-weak 2-s, 
then (GQP holds foripeC (Q; L™ r (Y x Z)) . See [9, Proposition 10] for the 
proof. 
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Instead of repeating here the main results underlying two-scale conver- 
gence, we find it more convenient to draw the reader's attention to a few 
references, see, e.g., [I], [7], [9] and [20] . 

However, we recall below two fundamental results. First of all, let 

y (0, T) = {ve L 2 (0, T; Hi (O; E)) : v' G L 2 (0, T; ff" 1 (O; E)) } . 
J 7 (0, T) is provided with the norm 

( 2 1 1 / 1 1 2 \ ^ 

IMIyfo.T) = W v h^(n.T:Hi(n^ + r \\T2(c)_T:H-i(n\) ( u G ^ (0,T)) 



l^(0,T) — \" v "L 2 (0,T;H](fy) IT 1 1 L 2 (0,T;H- 1 (H)) 

which makes it a Hilbert space. 

Theorem 3.1. Assume that 1 < p < oo and further E is a fundamental 
sequence. Let a sequence (u £ ) eeE be bounded in LP (Q). Then, a subsequence 
E' can be extracted from E such that (u £ ) £( - E , weakly two-scale converges in 
LP(Q). 

Theorem 3.2. Let E be a fundamental sequence. Suppose a sequence 
{ u e) £& E * s bounded in y(0,T). Then, a subsequence E' can be extracted 
from E such that, as E' 3 e — > ; 

u £ — > uo in y (0, T) -weak, 

u e — > uq in L 2 (Q) -weak 2-s, 

+ inL 2 (Q) -weak 2-s (1 < j < N) , 



du £ 8uq i du\ . r o 



dxj dxj dy 
where u G y (0, T), u\ G L 2 (Q; L 2 per (z; LL\ (Y) 

The proof of Theorem 13 . 1 1 can be found in, e.g., [I], [7], whereas Theorem 
2] has its proof in, e.g., [9] and [13] . 



3.2. A global homogenization theorem. Before we can establish a so- 
called global homogenization theorem for (|1.3p - (]1.6p . we require a few basic 
notation and results. To begin, let 

Vy = U G C£. (Y; R) N : (y) dy = 0, div y ij> =o| , 

y y = jw G H\ (Y; R) N : div y w =o} , 

where: C™ r {Y;R) = C°°(R N ;R) nC per (Y), div y denotes the divergence 

operator in Ry ■ We provide Vy with the (Y) N -norm, which makes it 
a Hilbert space. There is no difficulty in verifying that Vy is dense in Vy 
(proceed as in [144 Proposition 3.2]). With this in mind, set 

Fj = L 2 (0, T; V) x L 2 (Q; ( Z; . 

This is a Hilbert space with norm 



F 1 

IJT 



l V 0lli2(0,T;U) + l! V lllL2(Q;L2 er (Z;VV))) ' ' V = ( V 0> V l) G Fj 
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On the other hand, put 

F§° = V (0, T; V) x [V (Q; E) ® [C^ r (Z; R) ® Vy] ] , 

where C™ r (Z;R) = C°° (R; R) n C per (Z), C££. (Z;R) Vy stands for the 
space of vector functions w on My x IR T of the form 

w(y,r)= Xi(r)Vi(y) (r G R, y G R") 

with Xi G C~ r (Z; R), Vi G Vy , and where D (Q; R) (C™ r (Z; R) (g> Vy) is 
the space of vector functions on Q x R^ x R of the form 

i/>(x,t,y,r) = <Pi(x,t)wi(y,r) ({x,t)eQ, (y,r) € R N xR) 

finite 

with ^g23 (Q;R), Wj G C~ r (Z;R) <g> Vy. Since V is dense in V (see [IB| 
p. 18]), it is clear that Fq* 3 is dense in Fq. 
Now, let 

U = VxL 2 (n;L 2 per (Z;Vy)). 
Provided with the norm 

Hlu= (llvof+llvill^^.^^.^))) 2 (v= (v , Vl ) GU), 
U is a Hilbert space. Let us set 

ij^i •/ ■/ JnxYxZ \oxj dyjj \dxi dyi J 

for u = (uo, Ui) and v = (vq, vi) in U. This defines a symmetric continuous 
bilinear form on U x U. Furthermore, is U-elliptic. Specifically, 

(3.2) on(u,u) > a ||u||u (uGU) 

as is easily checked by using (jl.2p and the fact that j Y (x, y, r) dy = 0. 
Here is one fundamental lemma. 

Lemma 3.1. Under the hypotheses il.l\) - (T!?\) . The variational problem 
(3.3) 

' u G W (0, T) with u (0) = 0; 
u= (u ,ui) GFj : 
Jo (u' (t) , vq (i)) (it + J T ac (u (t) , v (i)) eft = J T (f (*) , v (t)) dt 
for all v = (v , vi) G F^ 

has at most one solution. 

Proof. Let v* = (vo, vi) G U and ip G V (]0,T[). By taking v =(p <8> v* in 
(13. 3p . we arrive at 

(3.4) (u / (t),vo)+an(u(i),v*) = (f(t),v ) (v* G U) 
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for almost all t € (0, T). Suppose that u* and u** are two solutions of (|3.3p 
with u* = (u*o,u*i) and u** = (u**o,u**i). Let u = u* — u** = (uo,Ui) 
with uo = u*o — u**o and u i = u *i — u**i- Let us show that u =0. By using 
(13.411 we see that u verifies: 



(3.5) (u' (i),v ) +a n (u(i),v*) = 
for all v* G U and for almost all t £ (0, T). But, by virtue of [I8j p. 261] 

||u (t)| 2 = 2(u' (t),u (t)) 
for almost all t 6 (0, T). Then, taking v t = u (t) in (|3.5p . we obtain by 

(3.6) j t \u (t)\ 2 + 2a\\u(t)\\l <0 

for almost all t £ (0,T). Integrating on [0,t] (0 < t < T), we get 

l u o(^)| 2 < for all i £ [0,7"] and ||u||pi < 0, thus u = and the lemma 
follows. □ 

We are now able to prove the desired theorem. Throughout the remainder 
of the present section, it is assumed that aij is Y-periodic for any 1 < i, j < 
N. 

Theorem 3.3. Suppose that the hypotheses of Lemma \2.1\ are satisfied. For 
< e < 1, let u £ be defined by tl.3\) - [T75\) . Then, as £ we have 

(3.7) u e — >• uo in W (0, T) -weak, 

(3.8) M M + M in L 2 (Q) _ weak 2 _ s (1 < jk < N) 
oxj oxj oyj 

where u = (uo,ui) (with uq = (uq) and u\ = (ufj) is the unique solution 

of U3i- 

Proof. By Proposition l2.2| we see that the sequences (Pe)o< e <i an d ( u e)o<e<i 
{u\, u ^) 0<£< i are bounded respectively in L? (Q) and W (0, T). Fur- 
ther, it follows from (|2.5p and (|2.6p that for 1 < A; < N, the sequence 
( n e)o<e<i * s bounded in y (0,T). Let £ be a fundamental sequence. Then, 
by Theorems I3.lti3.2l and the fact that W (0, T) is compactly embedded in 
L? (Q) , there exist a subsequence -E' extracted from I? and functions uo = 
K)i< fc <7v e W(0,T), Ul = K) 1<fc<JV € L 2 (Q;L 2 er (Z;H^ (Y;R) N 



and j? G L 2 (Q; L 2 er (Y x Z; R)) such that as^'se^O, we have ff377|> - ff378|) 
and 

(3.9) u £ — > uo in L 2 (Q)^ -strong, 

(3.10) p £ ->• p in L 2 (Q) -weak 2-s. 

But, by virtue of Lemma 13. 1\ the theorem will be entirely proved if we 
show that u = (uq,ui) verifies (|3.3p . In fact, according to (jl.4p . we have 
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div\iQ = and div y ux = 0. Therefore u = (uo, ui) G Fq. Let us recall that 
uo can be considered as a continuous function of [0, T] into H since W (0, T) 
is continuously embedded in C ([0, T] ; H). Let us show that uo (0) = 0. For 
v G V and tp G C 1 ([0, T]) with (/? (T) = and (/? (0) = 1, we have by an 
integration by part 

/ (u' e (t),v)<p(t)dt+ f (u £ (t),v)cp'(t)dt = -(u £ (0),v). 
Jo Jo 

According to (|1.6p , we have by passing to the limit in the preceding equality 
as E' 3 e -> 

/ (u (f),v) p(i)dt + / {u {t),v)(p'(t)dt = 0. 
Jo Jo 

Hence (uo (0) , v) = for all v G V, and as V is dense in H we conclude 
that uo (0) = 0. Now, let us check that u = (uo,Ui) verifies the variational 
equation of ()3.3f) . For < e < 1, let 



(3.11) 



■00 + eip\ with tp G V (Q; ] 



and 



01 GP(Q; 



oo (nr. ■ 

per \ i ■ 



Vy] , 



i.e., <& e (a?,i) = tf) (x,t) + etp 1 (x,t, |,|) for (sc,i) G Q. We have 3> £ G 
P(Q;R) Ar . Thus, multiplying (TO]) by * £ yields 



(3.12) 



/ T « (t) , <I> £ (t)) + J V (u £ (t) , <I> £ (i)) dt 
- f Q p £ div& £ dxdt = f Q (f (i) , $ £ (t)) dt. 



Let us note at once that 

i-T 



rp TV 

/ {u' £ (t),*e(t))dt = -J2 

Jo l=1 Jq 



dt 



+ e 



dxdt. 



Then by virtue of 
(3.13) " 

rT 



we have 

JV 



/ {u' £ (t),* £ (t))dt^-j2 [ 

Jo l=l Jq 

as E' 3 e — > 0. In fact, on one hand 



ui^-dxdt 
u dt 



N 

E 

z=i 



dt 



+ £ 



gv4 
if 



+ 



N 

E 



0" 



at 



dxdt + 



as £" 3 e — > 0, on the other hand 




xYxZ 



8r) 

dr 



dxdt 



dxdtdydr 
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f f f Q xYxZ u olfr dxdtd y dT = Iq u o (JJyxZ ^F d y dT ) dxdt = by virtue 
of the Y x Z-periodicity. The next point is to pass to the limit in (|3,12p as 
E' 3 e — > 0. To this end, we note that as E' 3 e — > 0, 



L 



f a £ (u E (t),& £ (t))dt^ [ a n (u(t),$(t))dt, 
o J o 

where $ = (V'ojV'i) (proceed as in the proof of the analogous result in 
[131 P-179]). Now, based on ()3.10f) . there is no difficulty in showing that as 
E' 3 e -> 0, 

/ p £ div<& £ dxdt pdivrp dxdtdydT. 

Jq J J JQxYxZ 

On the other hand, let us check that as e — > 

(3.14) F (f (t) , * £ (*)) dt -> / T (f (t) , (*)) 

Jo Jo 

Indeed, iff 6 L 2 (o, T; L 2 R) N ) (^T4"l) is immediate by using the classical 



fact that * £ ->■ i/; in L 2 (Q^-weak and §f? -> ^ in ^2 (g)^. weak 

(1 < j < iV) as e — ?> 0. In the general case, (|3.14p follows by the density of 

I? (o, T; L 2 (ft; R) N } in L 2 (o, T; H~ l (Q; R) N } . 

Having made this point, we can pass to the limit in (|3.12p when E' 3 e — > 
0, and the result is that 

(3 15) io T (u' (t) , (*)) dt + f* a n (u (t) , * (t)) eft 

- / Q p divjp dxdt = jf (f (f) , ^ (t)) dt, 

where po denotes the mean of p, i.e., po € L 2 (0, T; L 2 (Q; R)J and po (^j = 

/ fYxzP( x > t >y> T ) d y dT a - e - in G Q-> and wnere * = {^o^i)' V'o 

ranging over £> (Q; R)^ and ^ ranging over V (Q; R) £§ [C££ r (Z; R) <g) Vy] . 
Taking in particular i/> in D(0, T; V) and using the density of Fq° in Fj, 
one quickly arrives at (|3.3|) . The unicity of u = (uo, Ui) follows by Lemma 
I3.1[ Consequently, (|3.7p and (|3.8|) still hold when E 3 e — > 0. Hence when 
< e — > 0, by virtue of the arbitrariness of E. The theorem is proved. □ 

Now, we wish to give a simple representation of the vector function ui 
in Theorem I3.3I for further uses. For this purpose we introduce the bilinear 
form a on L 2 er (Z;Vy) X L 2 er (Z;Vy) defined by 

f f du k dv k 
a(u,v)= £ / / ^ ——dydr 



i,j,k=l 



YxZ vyj 
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for u = (u fc ) and v = (w fc ) G Lp er (Z; Fy). Next, for each pair of indices 
1 < i, k < N, we consider the variational problem 

( Xik^L 2 per (Z;V Y ): 

( 3 - 16 ) \ a (Xifc. w) = Ezli Jy x z a li^ d V dT 

{ for all w = (w?) G (Z; W) , 

which determines Xik m a unique manner. 

Lemma 3.2. Under the hypotheses and notation of Theorem \3. 31 we /iat>e 
(3-17) ui (x, t,y,r) = - ^ |j (x, t) Xik (y, r) 

i,k=l 1 

almost everywhere in (x,t,y,r) G Q x Y x Z . 

Proof. In (|3.3p . we choose the test functions v = (vo,vi) such that vo = 
and vi (x, t, y,r) = ip (x, t) w (y, r) for (x, t,y,r)£QxYxZ, where if G 

V (Q; M) and w G Lp er (Z; Vy). Then for almost every (x, t) in Q, we have 

(3.18) ( S ( Ul ' w ) = " £jj,*=l (fi, t) J J YxZ aij^dydr 
\ for all w G L 2 per (Z; Vy) . 

But it is clear that ui (x, t) (for fixed (x, t) G Q) is the unique function in 
Lp er (Z;Vy) solving the variational equation (|3.18|) . On the other hand, it 

is an easy exercise to verify that z (x, t) = — X)ffc=i t) Xik solves also 

(|3.18p . Hence the lemma follows immediately. □ 

3.3. Macroscopic homogenized equations. Our aim here is to derive a 
well-posed initial boundary value problem for (uo,po)- To begin, for 1 < 
hj, k,h < N, let 

Qijkh = hh J aij (y)dy ~YlJ J a ' il ^ ~ a ~ ( y > r ) dydT > 

where: Sp-h is the Kronecker symbol, Xjh = (x'jhj 1S defined by (|3.16p . To 
the coefficients qijkh we associate the differential operator Q on Q mapping 

V (Q) N into V {QY (V (Q) being the usual space of complex distributions 
on Q) as 

(3.19) 

N f) 2 h 

(Qz) k = - ]T q^khg^ (1 < k < N) for z = (z h ) ,z h GV' (Q) . 
i,j,h=l 1 J 

Q is the so-called homogenized operator associated to P e (0 < e < 1). 
Now, let us consider the initial boundary value problem 

(3.20) ^£ + Qu + gradp = f in Q = Ox]0,T[, 

at 
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(3.21) 
(3.22) 



divuQ = in Q, 
u = on dfix]0,T[, 



(3.23) u (0) = in Q. 

Lemma 3.3. The initial boundary value problem t3.20\ )- [3.23\) admits at 
most one weak solution (uo,po) with 

u G W (0, T) and p G L 2 (0, T; I? (U; R) /R) . 

Proo/. If (u ,po) G W (0, T) x L 2 (0, T; L 2 (fi; R)) verifies (IH^UD - (13^31 . then 
we have 

/cT K (*) > v o (*)) dt + EZ,k,h=i Iq Hjkh^^dxdt 
= / T (f(t),v (t))dt 

for all v G L 2 (0, T; V). From the previous equality, one quickly arrives at 
(3.24) 

Jo K (t) , v (t)) + Efi,fc=i / / Jgxyxz a ii (Sj + fj") ^dxdtdydr 

where «J (x, t, y, r) = - £f fc=1 §J ( x ' *) (2/' r ) for ( x > *>!/>t) G QxFxZ. 

Let us check that u = (u , ui) (with ui (x, t, y, t) = - 22 i)k=1 (x, t) Xik \V, T ) 
for (x,t,y,r) & Q x Y x Z) satisfies (|3.3p . Indeed, we have 

for all vi = (i^) G L 2 (Q; L 2 er (Z; Vy)) , since ui(x,i) verifies f|3. 18|) for 
(x,t) G Q. Thus, by (I3"^4l - (l3%| ). we see that u = (u ,ui) verifies (f3T3]) . 
Hence, the unicity in (|3.20p - (|3.23p follows by Lemme I3TT1 □ 

This leads us to the following theorem. 

Theorem 3.4. Suppose that the hypotheses of Theorem are satisfied. 
For each < e < 1, let (u E ,p e ) G W (0, T) x L? (0, T; L 2 (fl;R) /R) 6e 
defined by $1. &)) - $!. b}) . Then, as e — > 0, we have u e — > uq in W(0,T)- 
weafc and p £ — > po in L 2 (0,T;L 2 (Q)) -weak, where the pair (uo,po) ^es in 
W (0, T) x 1? (0, T; L 2 (O; M) /R) and is tte umoue solution o/ M) - JI^) . 

Proof. Let -E 1 be a fundamental sequence. As in the proof of Theorem 13. 3| 
there exists a subsequence I?' extracted from E such that as E' 3 e — >• 0, we 
have (|37f l> -(|33 ]l and ([3T0]) with u = (u ,ui) G Fj and u (0). Then, from 
(|3.10p we have p £ — > po in L 2 (0, T; L 2 (Jl))-weak when E' 9 e — > 0, where po 
is the mean of p. Hence, it follows that p G L 2 (0, T; L 2 (Q; R) /R) . Furher, 

PUD holds for all * = (?/> , Vi) G P (Q; ffxP (Q; R)® [C~ r (Z; R) ® Vy] . 
Then, substituting (|3.17p in (|3. 15[) and choosing therein the <&'s such that 
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ij>i = 0, a simple computation leads to (|3.20p with evidently (|3.21|) - (j3.23p . 
Hence the Theorem follows by Lemma 13.31 since E is arbitrarily chosen. □ 

Remark 3.2. The operator Q is elliptic, i.e., there is some ao > such 
that 

N N 

Qijkh^tkCjh > «o ^ \ikh? 

i,j,k,h=l k,h=l 

for all £ = with ^ 6 M. Indeed, by following a classical line of argu- 
ment (see, e.g., [2]J, we can give a suitable expression of qijkh, viz. 

Qijkh = a [Xik - ^ik, Xjh - 77 jh) , 
where, for each pair of indices 1 < i, k < N , the vector function tvh- = 

( 7r ijfe>--> 7r ife) : R y R is a iven b v Kkiv) = yrfkr (r= h-, N ) f° r v = 

(yi, 2/jv) G • Hence, the above ellipticity property follows in a classical 
fashion. 
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